Ab initio calculation of the lattice distortions induced by substitutional Ag 

impurities in alkali halide crystals. 
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An ab initio study of the doping of alkali halide crystals (AX: A = Li, Na, K, Rb; X = F, CI, 
Br, I) by ns^ anions (Ag~ and Cu~) is presented. Large active clusters with 179 ions embedded in 
the surrounding crystalline lattice are considered in order to describe properly the lattice relaxation 
induced by the introduction of substitutional impurities. In all the cases considered, the lattice 
distortions imply the concerted movement of several shells of neighbors. The shell displacements are 
smaller for the smaller anion Cu~, as expected. The study of the family of rock-salt alkali halides 
(excepting CsF) allows us to extract trends that might be useful at a predictive level in the study 
of other impurity systems. Those trends are presented and discussed in terms of simple geometric 
arguments. 
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I. INTRODUCTION 

Most of the luminescent materials presently used in 
several technological applicationsB involve the doping of a 
pure ionic crystal, that is substitution of some of the ions 
by other ions with specific absorption-emission character- 
istics. The fine details of the absorption-emission spectra, 
as well as the efficiency and resolution of the scintillator, 
are determined by the system-specific embedding poten- 
tial acting on the impurity, which is in turn sensitive 
to the distortion induced by the impurity on the crystal 
lattice. Thus, a theoretical understanding and accurate 
determination of those distortions is of paramount im- 
portance, moreover if we realize, that their experimental 
measurement is a difficult taskjjlj 

Two main methods are applied nowadays to model 
impurity systems: supercell techniques, that exploit the 
convenience of the Bloch theorem by periodically dup 
eating a finite region of the crystal around the impurit; 
and the cluster approach, in which the doped crystal is 
modeled by a finite cluster centered on the impurity and 
embedded in a field representing the rest of the host 
lattice. This cluster approach is the one chosen in the 
present study, and has been used in the past to studyp 
geometrical and optical properties of doped crystals 
The cluster (active space) can be studied by using stan- 
dard quantum- mechanical methods. The rest of the crys- 
tal (environment) can be described in several ways. In 
the simplest and most frequently used approach, the en- 
vironment is simulated by placing point charges on the 
lattice sites, but this procedure has to be improved in 
order to obtain a realistic desc.rijpt jnUpXLLihe Jaitice dis- 
tortions around the impurityliayil3'BNEa'E3y Model 
potentials have been developed to represent the effects of 
the environment on the active cluster, that include at- 
tractive and repulsive quantum-rnechanical terms aside 
from the classical Madelung termj^ but a problem still 



remains: the large computational cost of conventional 
molecular orbital (MO) calculations prevents from per- 
forming an exhaustive geometrical relaxation of the lat- 
tice around |4lie,riHmT«i|tyi. In the most accurate MO 
calculations ,ll3E3lHa'EZl'EjH only the positions of the ions 
in the first shell around the impurity are allowed to re- 
lax. However, geometrical relaxations far beyond the 
first shell of neighbors can be expectejd, r-Ln fact, re- 
cent semiempirical simulations of solids ir.Trp performed 
employing phenomenological potentials have shown 
the importance of considering appropriate large-scale lat- 
tice relaxations in the study of a variety of intrinsic and 
extrinsic defects in ionic crystals. As we will show below, 
the systems under study in this paper can not be prop- 
erly described by simply considering an expansion of the 
first shell of neighbors around the impurity. 

In this contribution we report theoretical calculations 
of the lattice distortions induced by Ag^ and Cu^ sub- 
stitutional impurities in 16 alkali halide crystals with the 
rock-salt structure, namely all those noncontaining ce- 
sium. For this purpose we use the ab initio Perturbed Ion 
(PI) model,E2rE^ which circunvents the problems men- 
tioned above: (a) The active cluster is embedded in an 
environment represented by the ab initio model poten- 
tials of Huzinaga et alB3; (b) The computational simplic- 
ity of the PI model allows for the geometrical relasatiM, 
of several coordination shells around the impurity.Mrl'Ell 
Moreover, it allows us to study a whole family of systems 
in order to look for systematic trends that might be use- 
ful in later theoretical studies of doped crystals similar 
to those here considered. 

The remainder of this paper is organized as follows: 
In Section we describe the active cluster which has 
been used to model the doped systems. In Section [II we 



present and discuss the results of the calculations, and 
Section summarizes the main conclusions. 



1 



II. CLUSTER MODEL 

The ab initio Perturbed Ion model is a particular ap- 
plication of the the ory of electronic separability of Huzi- 
naga and coworkera^Ej to ionic solids, in which the basic 
building blocks are reduced to single ious.. The PI model 
was first developed for perfect crystals .CJ Its application 
to the study of im purity centers in ionic crystals has been 
described in refs. 15|,|l8|,^ and we refer to those papers 
for a full account of the method. In brief, an active cluster 
containing the impurity is considered, and the Hartree- 
Fock-Roothaan (HFR) equationst3 for each ion in the ac- 
tive cluster are solved in the field of the other ions. The 
Fock operator includes, apart from the usual intra-atomic 
terms, an accurate quantum-mechanical crystal poten- 
tial and a lattice projection operator which accounts for 
the energy contribution due-i;o the overlap between the 
wave functions of the ionsH The atomic-like HFR so- 
lutions are used to describe the ions in the active clus- 
ter in an iterative stepwise procedure. The wave func- 
tions of the lattice ions external to the active cluster are 
taken from a PI calculation for the perfect crystal and 
are kept frozen during the embedded-cluster calculation. 
Those wave functions are explicitely considered for ions 
up to a distance d from the center of the active cluster 
such that the quantal contribution from the most distant 
frozen shell to the effective cluster energy is less than 
10~^ hartree. Ions at distances beyond d contribute to 
the effective energy of the active cluster just through the 
long-range Madelung interaction, so they are represented 
by point charges. At the end of the calculation, the ionic 
wave functions are selfconsistent within the active cluster 
and consistent with the frozen description of the rest of 
the lattice. The intraatomic Coulomb correlation, which 
is neglected at the Hartrce-Fock level, is computed as 
a correction by usiag_Lhe Coulomb-Hartree-Fock (CHF) 
model of Clementi.Ht; 

In a previous workE^I we employed several active clus- 
ters of increasing size and with different embedding 
schemes to describe the scintillator system Tl+iNal. That 
study was undertaken in order to find the necessary re- 
quirements that a cluster model has to fulfill in order 
to describe properly a doped crystal. Here we just de- 
scribe the best cluster model between those studied in 
ref. This active cluster, shown in Figure 1, has 179 
ions which correspond to the central impurity (Ag~ or 
Cu~) plus twelve coordination shells. Those ions are fur- 
ther split up into two subsets. One is formed by the 
central impurity plus the first four coordination shells, 
having a total of 33 ions, and both the wave functions 
and positions of the ions in this subset are allowed to 
relax. The lattice positions of the other 146 ions of the 
active cluster are held fixed during the calculations but 
their wave functions have been selfconsistently optimized. 
This is done so that the connection between the region 
where distortions are relevant and the rest of the cryja^ 
tal is as smooth as possible. In our previous studyEll 



we showed how an unphysically abrupt connection be- 
tween those two regions fails in describing properly the 
lattice distortions induced by the impurity. The ions in 
the interface region can respond to those distortions by 
selfconsistently adapting their wave functions to the new 
potential, and thus contribute to build a more realistic 
(selfconsistent) environment. The geometrical relaxation 
around the impurity has been performed by allowing for 
the independent breathing displacements of each shell 
of ions, and minimizing the total energy with respect to 
those displacements until the effective cluster energies sf«. 
converged up to 1 meV. A downhill simplex algorithmE^I 
was used. For the ions we have used larse, STO basis 
sets, all taken from Clementi-Roetti tablesH 

The cluster used in this work has been shown to be-self- 
embedding consistent for Nal in our previous work.LJ By 
this we mean that if the pure crystal is represented by this 
cluster model (that is, if the central impurity is replaced 
by the halogen ion corresponding to the pure crystal), 
the results of the cluster model calculations closely agree 
with those from a PI calculation for the pure crystal, 
where all cations (or anions) are equivalent by transla- 
tional symmetry. The same is true for all the family of 
alkali halide crystals considered here. Nevertheless, the 
self-embedding consistency is never complete. In order to 
supress systematic errors from the distortions calculated 
with the cluster method, the radial displacements of each 
shell have been calculated using the following formula: 



AR, = R,{Imp- : AX) - Ri{X- : AX), 



(1) 



where R^ (i=l, 2, 3, 4) refer to the radii of the first, 
second, third, and fourth shells around the impurity in 
the AX crystal, A = Li, Na, K, Rb, X = F, CI, Br, I, 
and Imp" = Ag~, Cu~. Thus both systems (pure and 
doped crystals) are treated in eq. (1) on equal foot with 
the cluster model, and not with different methodologies, 
and the calculated distortions are free from that poten- 
tial source of error. Also, in order to have the correct 
Madelung potential at the impurity site, the calculations 
have been performed by employing the experimental lat- 
tice constantg23 to describe the geometrically frozen part 
of the crystals. 

The only terms omitted in our description are the dis- 
persion terms (coming from interatomic correlation) and 
relativistic effects for the heavy ions. Although the im- 
portance of both effects increase with atomic number, 
they are not crucial for the structural properties of the 
systems studied here. Specifically, Martin Pendas et a/.E3 
have shown that the PI method gives lattice constants 
and bulk moduli in close agreement with experimental 
results for all alkali halides. The properties of these crys- 
tals under the influence of an applied external pressure, a 
situation where the importance of interatomic correlation 
effects increases, are also properly reproduced. 



2 



III. RESULTS AND DISCUSSION 

The calculated distortions, collected in Table I, are the 
main quantitative result from our study. For visualiza- 
tion of the trends, however, it is better to display the 
results in a figure, and this is done in Fig. 2, where we 
have plotted the distortion of each of the four shells in 
terms of the empirical cationic radii extracted from ref. 
|62| . Those points corresponding to the same anion have 
been joined with a line to guide the eye. The figure con- 
tains only the results for Ag~ because the trends are the 
same in the case of the Cu^ impurity. Next we describe, 
shell by shell, the general trends in Fig. 2: 

First shell. This shell is formed by 6 cations in (^,0,0) 
crystallographic sites, and undergoes an expansion, as 
might be expected from the larger size of Ag~ compared 
to the halogen anions. The impurity anion pushes the 
neighbor cations to make room for itself in the lattice. 
The expansion is substantial, with percentage values be- 
tween 9 and 15 %. In the F salts that expansion is larger 
the larger the cation size, but this trend is violated in 
the CI, Br, and I crystals. If we fix the cation, for K 
and Rb salts the expansion is larger the smaller the an- 
ion (notice that the anion size increases in the order F~ , 
Cl^, Br^, I^). This rule is inverted in the case of Li 
salts, whereas Na salts constitute an intermediate case. 
It should be recognised, nevertheless, that the expansion 
is almost independent of the halogen element in the Li 
and Na salts. 

Second shell. The displacement of the second shell, 
formed by 12 anions at (^,^,0) positions, is always a 
small contraction. If we fix the anion, the contraction 
is larger (absolute value) the larger is cation size. If the 
cation is fixed, for Na, K and Rb salts the contraction 
is larger the smaller the anion size. Again this trend is 
inverted for Li salts. 

Third shell. This shell, formed by 8 cations at (^,^,5) 
sites, experiences a small expansion. If we fix the anion, 
the expansion is smaller the larger the cation size. If the 
cation is fixed, there is not a definite trend. In the case 
of Li and Na salts the expansion increases with the anion 
size (Lil is an exception) . This trend is partially inverted 
in the case of Rb salts, and for K salts all the expansions 
are almost identical. 

Fourth sell. This shell, formed by 6 anions at (1,0,0) 
positions, experiences an expansion. That expansion in- 
creases with cation size if we fix the anion. If the cation 
is fixed, in K and Rb salts the expansion is smaller the 
larger the anion size. In Na salts, Nal is again an excep- 
tion to this general rule, whereas in Li salts the expansion 
is almost constant. 

In the following we try to find some rationalization 
for the calculated trends. The working rule still in use 
nowadays stating that the expansion of the first coordi- 
nation shell can be approximated by the difference be- 
tween the ionic radii of the impurity and the substituted 
ion is somewhat misleading. First of all, the ions in a 



crystal are not hard spheres, but weakly overlapping soft 
spheres, so one cannot use values for the ionic radii of 
ions in vacuum in order to predict lattice distortions ac- 
curately. In particular, the size of an anion may vary 
in a nonnegligible way from crystal to crystal. To in- 
vestigate this, we show (< r^ >)^/^ for F^ in fluoride 
crystals in Table II, where the expectation value is taken 
over the outermost orbital of the anion (2p), and is cal- 
culated from the crystal-consistent ionic wave functions 
obtained through a PI calculation on the pure crystals. 
Tx — (< r^ >)^/^ can be taken as a rough measure of 
the anion size. We also show the analogous quantity for 
the 5s orbital of the Ag~ impurity in fluorides. The size 
of the F^ anion varies by a maximum of 2 %, small com- 
pared with the size variation of the Ag" anion (7 %). 
Ag~ is more compressible than the halogen anions be- 
cause its outer electronic shell is an s-shell, while it is 
a p-shell for the halogens. The same can be said of all 
the halogens. As the size of the cation decreases, Ag~ is 
more compressed by the crystal environment. This shows 
that the standard ionic radii cannot be used for predict- 
ing distortions, because their values are a genuine output 
of the selfconsistent process. Nevertheless, though they 
are not useful for accurate predictions, physical insight 
tells us that the distortions should be correlated with the 
size of the ions, in the sense that one always expect that 
larger impurities induce larger distortions. These should 
be useful at least at a qualitative level. 
In Table III we show the differences 

S = riAg- ■.AX)^rx{X- -.AX), (2) 

where r(Ag~ : AX) is the radius of the Ag^ impurity in the 
AX crystal, and rj(:(X~:AX) is the radius of the halogen 
anion. The differences are expected to be related with 
the flrst-shell expansions Ai?i. The values of 6 do not 
show quantitative agreement with those of Ai?i . If the 
anion is fixed, S increases with the cation size, which is 
consistent with the main trend in Ai?i of fluorides, but 
does not explain the behavior of other halides. The main 
trends discussed when the cation is fixed are reproduced 
for the K and Rb crystals but not for the others. From 
Figure 1 we can see that the exceptional crystals, con- 
cerning the trends in Ai?i, are LiCl, LiBr, Lil, and Nal. 
The peculiar feature of those four cases is that anions 
are much larger than cations: the ratios ta/^x are the 
smallest in the family of alkali halides. In_a study of the 
structures of small alkali halide clusters)^ it was found 
that materials with small r^^/rx have a cluster growing 
pattern different from the rest. Specifically, those sys- 
tems showed a marked tendency towards ring-like struc- 
tures, whereas the others adopt fragments of the rocksalt 
lattice as their minimum energy structures. The main 
reason is that anion-anion repulsions are much more im- 
portant when the ratio ta^ix is small. In a recent study, 
A. Martin Pcndas et aZ.Ej havej-applied the atoms-in- 
molecules (AIM) theory of BadeiEj to study the topology 
of the electron density in crystals. They found that the 
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whole group of alkaH hahdes with rocksalt structure can 
be divided up into three topological families, called Ri, 
Bi, and B2 in that paper, and that the ionic radii are 
the topological organizers. The Ri family contains KCl, 
NaF, KF, RbF, RbBr, RbCl and all the cesium halides. 
The Bi family contains KI, KBr, LiF, Rbl, NaBr, NaCl 
and NaT, and the B2 family contains LiCl, LiBr and Lil. 
There are constant ia/'^x lines that isolate each fam- 
ily. The largest values of r^/rx are found for family Ri 
and the smallest for family B2, with intermediate val- 
ues for Bi. NaT is so near the B2 region that it is not 
surprising that it behaves in Fig. 2 like the elements 
of the B2 family. In the AIM theory the critical points 
of the electron density scalar field are classi: 
clei, bond points, ring points and cage points, 
bond point is found between two nuclei, a bond is estab- 
lished between the corresponding atoros,. In the Ri fam- 
ily there are just anion-cation bonds.E3 We have found 
that for the (undoped) crystals of this family, the anion- 
anion overlap is at least one order of magnitude smaller 
than the anion-cation overlap. In the crystals of the Bi 
and B2 families there are bond critical points between 
anions, so that the effective local coordination is 6 for 
cations and ISJor anions (6 anion-cation and 12 anion- 
anion bonds) .l3 In the Bi crystals (excepting Nal) we 
have found that the anion-anion overlap is smaller but of 
the same order of magnitude than cation-anion overlap. 
In the B2 crystals, and also in Nal, anion-anion overlap 
is the largest contribution to the repulsive interactions, 
and thus anion-cation contacts are less important. The 
cations of the B2 family occupy the interstitial holes left 
in the anionic fee sublattice. It is then not surprising 
that when the cation-anion overlap is not so important, 
the expansion of the cation shell is an exception to the 
general trends; it is, in fact, nearly constant for the B2 
family. 

Let us turn to discuss the distortion of the second and 
third shells. In all cases, the second shell suffers a con- 
traction of a small magnitude compared to the large ex- 
pansion of the first shell. In Figure 1 one can see that 
the radial outward motion of the first cation shell is not 
going to affect much the positions of the twelve anions 
of shell 2, so it becomes understandable that the anions 
of that shell move little. The small contraction of the 
second shell optimizes the Madelung energy around the 
impurity and also serves to pack more efficiently the ions 
in response to the outward motion of the cations. The 
quantitative trend of that contraction is understood with 
reference to the three topological families discussed in 
the previous paragraph and their relation to the anion- 
anion overlap: the contraction is largest for those crys- 
tals where the anion- anion overlap is small (Ri crystals), 
intermediate when that overlap begins to count (Bi crys- 
tals), and finally, it is lowest for the B2 crystals, where 
anion- anion contacts are important. This explains the 
trends observed: if the anion is fixed, the contraction 
is larger the larger the cation size, because the Ag~-X~ 
overlap decreases with increasing cation size. On the 



other hand, if the cation is fixed, in Rb, K, and Na salts 
the contraction decreases with increasing anion size, be- 
cause the Ag^-X^ overlap increases with anion size. But 
in Li salts the trend is inverted because Ag~-X~ over- 
laps decrease with anion size. The distortion of the third 
shell is not directly related to the introduction of the im- 
purity, as the overlap between the cations of that shell 
and Ag~ is very small (see Fig. 1). The expansion of 
this shell appears to be again of a purely electrostatic 
origin. The relative values of the displacements (R3 - 
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are very small, less than 1 % except in 
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LiF. 

The Ai?4 displacements, always an expansion, proceed 
along the same crystallographic direction as the Ai?i 
displacements. Thus, the expansion is clearly induced 
by the expansion of the first shell. The relative dis- 
placements (R4-R4'^'''**''')/R4'^^'**''' adopt values between 
1 % and 5%, compared to values of 9-15 % for (Ri- 
j^crystaiyj^crystai ^ Thcsc numbcrs indicate that the ex- 
pansion of the (ninth) shell formed by six cations at 
(|,0,0) is not expected to be higher than 1 %. In the 
K and Rb salts AR4 is larger when Ai?i is larger, while 
in Li and Na salts Ai?4 and ARi are both almost inde- 
pendent of the anion, so the displacements of the first 
and fourth shells are correlated. If the anion is fixed, the 
expansion increases with the cation size, and no special 
behaviour is observed in the cases of LiCl, LiBr, Lil and 
Nal. 

We conclude that the lattice relaxation around the sub- 
stitutional impurity in the alkali halides involves the con- 
certed movement of several coordination shells. However, 
it is not yet clear from the results presented up to this 
point whether the lattice relaxations of the second, third 
and fourth shells have a substantial infiuence on the en- 
ergy of formation of the defect. At low pressure and 
temperature conditions, the formation of the impurity 
centers should be discussed in terms of the internal en- 
ergy difference for the exchange reactionO 



{X- : AX)s + Imp- ^ {Imp- : AX), + X" 



(3) 



where the s and g subindexes refer to solid and gas 
phases, respectively. In order to establish the importance 
of the relaxation of the lattice beyond the first coordina- 
tion shell, we have calculated the formation energy AH of 
the defects by employing two different models for the ac- 
tive cluster, shown in Fig. 1: one of them is that formed 
by 179 ions described in section II; the other includes just 
four coordination shells around the central ion (a total of 
33 ions), and only the positions of the six ions in the first 
coordination shell are allowed to relax. The results are 
shown in Table IV. We see that according to the small 
cluster model the energy of formation of the defects is 
always positive, that is none of the impurity centers are 
stable centers. Enlarging the cluster size to include a self- 
consistent treatment of 179 ions and extending the lattice 
relaxation up to the fourth coordination shell induces a 
huge stabilization of all the impurity centers. The trend 
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of AH is simple. For a given impurity (Ag~ or Cu~) 
the heat of formation deereases by increasing the atomie 
number of the cation (alkah) or of the anion (halogen). 
The positive value of AH in the calculation for the small 
cluster is understandable: the impurity simply pushes its 
neighbor cations producing a high elastic strain energy 
since the rest of the lattice is not allowed to respond. In 
the second model three more shells are allowed to move 
in response to that initial stress and the relaxation lowers 
the elastic energy so much that the electronic contribu- 
tions turn AH negative in all cases except Ag~:LiF and 
Ag~:NaF (AH is nearly zero in Ag~:LiCl). It is useful 
to notice that ARi, the displacement of the first shell, 
is lower for the first cluster model compared to the sec- 
ond. This means that due to the constraints imposed by 
the first model the atoms of the first shell are unable to 
reach their preferred equilibrium positions in the presence 
of the impurity, a fact that is consistent with the large 
calculated AH. The change of sign in AH can be inter- 
preted as suggesting that most of the elastic relaxation 
of the lattice has been accounted for and that allowing 
for the elastic relaxation of more shells will have a minor 
effect. The two cases with a positive AH, Ag~:LiF and 
Ag~:NaF, are still intriguing. These two crystals have 
the smallest lattice parameters within the whole family 
studied here, and it is conceivable that the elastic effects 
will be largest. The question if the relaxation of more 
coordination shells is able to stabilize those two systems 
deserves further investigation. 

As indicated above the difference in the heats of forma- 
tion given by the two models, AH(model 2) - AH(model 
1), gives a measure of the lowering of elastic strain when 
more coordination shells are allowed to relax. Prom Ta- 
ble IV one can verify that, for a given crystal, this energy 
is essentially independent of the impurity, while both 
AH(model 1) and AH(model 2) depend on the impurity. 
This confirms our interpretation of the effect of allowing 
for the elastic relaxation of several shells: that relaxation 
is mainly a host effect. 

IV. SUMMARY 

We have reported a study of the local lattice distor- 
tions induced by substitutional Ag~ and Cu~ impurities 
in the family of alkali halide crystals excepting those con- 
taining cesium. For this purpose, the ab initio Perturbed 
Ion (PI) model has been used. A large active cluster of 
179 ions, embedded in an accurate quantum environment 
representing the rest of the crystal, has been studied. The 
local distortions obtained extend beyond the first shell of 
neighbors in all cases. Thus, the assumptions frequently 
employed in impurity calculations, which consider the ac- 
tive space as formed by the central impurity plus its first 
coordination shell only, should be taken with some care. 
Distortion trends have been identified and discussed. The 
first coordination shell (cations) around the impurity ex- 



periences an expansion as a consequence of the larger size 
of the impurity anion compared to the halogens. That 
expansion is larger for the Ag^ than for the Cu^ impu- 
rity, also because the first anion is larger than the second. 
The trends can be qualitatively explained by considering 
the difference in size between the impiirity and the sub- 
stituted anion in all cases except in those crystals with a 
very small size ratio between cation and anion. In those 
cases, that is for LiCl, LiBr, Lil and Nal, anion-anion 
contacts are important. Those four materials have been 
found to exhibit special behavior in a number of previous 
studies involving crystals and clusters. The contraction 
of the second shell as well as the expansion of the third 
shell are small and arise from a combination of electro- 
static and packing origins. The fourth shell experiences 
a substantial expansion as a consequence of the direct 
pushing induced by the expansion of the first shell. The 
analysis of the energies of formation of the defects clearly 
shows that elastic relaxation of several coordination shells 
around the impurity is necessary in the modeling of these 
materials since this affects even the sign of the energy of 
formation. 
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Captions of Tables 

Table I. Radial displacements ARi (see eq. (1)), in 
A, of the first four shells of ions around the silver and 
copper impurities. 

Table II. (< >)^^^, where the expectation values 
are taken over the outermost orbital of the F^ anion in 
pure alkali fluorides and of the silver anion in Ag~-doped 
alkali fluorides. All quantities in A. 

Table III. Difference of radii between Ag~ and the 
substituted anion (see eq. (2)), in A. 

Table IV. Formation energies (in eV) of copper and 
silver substitutional centers in different alkali halide host 
lattices, calculated employing two different models for 
the active cluster. First row: the active cluster contains 
33 ions, and only the positions of the ions in the first 
coordination shell are allowed to relax. Second row: the 
active cluster contains 179 ions, and the positions of the 
ions in the first four coordination shells are allowed to 
relax. 

Captions of Figures 

Figure 1. The active cluster (ImpA92X86)^''~ em- 
ployed to represent the region around the impurity, where 
Imp=Ag, Cu; light spheres are cations and dark spheres 
anions. The core of the cluster, formed by the four first 
coordination shells, which are allowed to breath, is also 
indicated separately. 
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Figure 2. Shell distortions AR^ (1=1,2,3,4) around 
Ag~, plotted as a function of cation size. 
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